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Abstract 

We investigate the possible existence of nonradiating motions of sys- 
tems of point charges, according to classical electrodynamics with retarded 
potentials. We prove that two point particles of arbitrary electric charges 
cannot move for an infinitely long time within a finite region of space 
without radiating electromagnetic energy. We show however with an ex- 
ample that nonradiating accelerated motions of systems of point charges 
do in general exist. 



1 Introduction 

It has already been proved since a long time (see [U [3 [31 |1] and the refer- 
ences given in j2] to the previous works of Herglotz, Sommerfeld, Hertz and 
Schott) that there exist nontrivial motions of extended electric charge distri- 
butions which do not radiate energy according to the classical theory of elec- 
tromagnetism. Some of these motions refer to rigid charged bodies, and were 
originally associated with classical extended models of elementary particles. In 
more recent times nonradiating sources have gained renewed popularity in con- 
nection with the study of the inverse problem in wave equations, namely, the 
problem of reconstructing a source when the radiation emitted or scattered by 
it is known (see for instance [SUHIIZ] and references therein). 

Most of the existing literature in this field deals with spatially extended 
sources with a monochromatic dependence on time. In the present paper we 
will instead look for nonradiating systems made of point charges, with no a 
priori restriction on their possible motion. Hence our mathematical approach 
to the problem of nonradiating sources will be completely different, and in some 
sense complementary to that usually followed. From a fundamental point of 
view, an obvious motivation for our study comes from the fact that, according 
to present-day standard theories of microscopic physics (with the exception 
of string theory), elementary electric charges in nature are actually point-like. 
Therefore, in these theories, continuous charge distributions can only serve as 
useful approximations for the description of macroscopic bodies. 

The results of the present investigation may be relevant in connection with 
the search for classical models of atomic systems. A classical atom is in fact 
usually described as an isolated system of moving point charges (nucleus and 
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electrons). Since the absence of radiation is a necessary condition for the sta- 
bility of the atomic ground state, the formulation of a classical atomic model 
can only be possible provided that nonradiating motions of point charges in- 
deed exist, and that they are compatible with suitable dynamical laws taking 
into account radiation reaction (one may adopt for instance the classical third 
order Lorentz-Dirac equation [SI [3] , or its approximated version of second order 
with respect to time [ini HI])- There is a widespread belief that all solutions of 
these dynamical equations do actually radiate, and that classical physics there- 
fore cannot account for the stability of atomic systems. Nevertheless, attempts 
to describe atomic physics by making only use of the laws of classical elec- 
trodynamics have still recently been undertaken [12j . Moreover, nonradiating 
motions have been found in the dipole approximation for infinite regular ar- 
rays (with arbitrary lattice parameter) of point-like charged oscillators obeying 
to the Lorentz-Dirac equation with retarded mutual electromagnetic interac- 
tions [I3t il4j . It would therefore be interesting, in our opinion, to establish in 
which cases the impossibility of nonradiating motions of a finite number of point 
charges can indeed be proved with rigorous mathematical arguments. 

In this work we shall not be directly concerned with the dynamical equations 
which determine the motion of the particles when the field acting on them is 
known. Our aim will be simply to study the restrictions which are imposed on 
any arbitrary motion of point charges by the condition of vanishing radiation. 
This radiation will be calculated by making use of the usual retarded electro- 
magnetic potentials which arise from Maxwell equations with point-like sources. 
Particles moving with constant velocities of course do not radiate electromag- 
netic energy, and their distance from any fixed point in space increases to infinity 
whenever their velocities are not zero. For the case of one or two particles with 
arbitrary charges, we will prove that there do not exist nonradiating motions 
which are bounded in space for all times. However, we will find that a pair of 
point particles of equal charges do not radiate while they move on a straight 
line in opposite directions, with their spatial coordinates varying as the square 
root of time. This is an indication of the fact that accelerated nonradiating 
motions of systems of point charges do in general exist. The search for other 
nontrivial examples involving a number of particles greater than two could be 
an interesting matter for further investigation. 

2 The condition of zero radiation 

Let us formulate our problem in more precise mathematical terms. We shall 
consider smooth motions of N point particles, with nonvanishing electric charges 
(7i , . . . , gjv of arbitrary magnitudes and signs. We denote by the three-dimensional 
vector Zi{t) the position of particle i at time t, with respect to a fixed cartesian 
coordinate system. We suppose that the motion of these charges is confined 
within a finite region of space. This means that there exists a fixed length 
i > such that 

\z,{t)\ < L for all t and alH = 1, . . . , iV . (1) 

The condition of no radiation is expressed as the requirement that the flux of the 
Poynting vector generated by the charges, calculated through any large spherical 
surface of radius L, vanishes at all times. Apart from the conditions just 
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mentioned, the motions considered will be a priori of the most general possible 
type. 

The retarded electric field generated at space-time point (x, t) by point par- 
ticle i is [T5] 



E,(x,t) 



i?,(l-n, -vO^ 



n, • V,; 



Ri 



n,; - V,: 



1 - rij • Vi 



(2) 



Here Ri — \x — Zi{ti)\ and = [x — Zi{ti)]/Ri (so that |ni| = 1), where ti is 
the retarded time of particle i, which is defined as an implicit function of (x, t) 
by the equation 

t-U = \jc-z,{U)\. (3) 

In ^ we put also v.; = Vi{ti) ~ dzi{ti) / dti and = Vi(ii) — dvi{ti)/dti. The 
magnetic field can then be expressed as 



Bi(x, i) = rij X Ej(x,t) . 



(4) 



Note that in our units the speed of light is 1. 

Let us evaluate the retarded field ([2]) at a point x such that |x| = R. In the 
limit i? ^ CXI we have = n -I- 0{R~^), where the unit vector n = x/i? defines 
the particular direction considered. We also have Ri = R — n- Zi{ti) + 0{R"^). 
Therefore, if we call t^, the time at which the field is evaluated at x, using Eq. 
Q and neglecting infinitesimal terms we can write 

= t + n • Zj(ti) , (5) 

where t ^ tji — R. It follows that for R ^ oo Eq. ^ can be simplified as 



E,(x,t«) 



i?(l-n-v,)2 



n ■ Vt 

1 — n • V, 



■(n- V,) 



+ 0{R-^), (6) 



where and are evaluated at the time U which is implicitly defined by Eq. 

It is convenient for our purposes that in this equation t be considered 
independent of R. This means that the retarded time U is also independent of 
R, whereas tji = t + R must increase with R for fixed t. 

For the total fields generated by the system of particles we have E = _R^^E + 
OiR-^), B = R-^n X E 0{R-^), where 



E 



N 



1 - n • V,; 



(n - V,) 



(7) 



is a quantity independent of R. The power radiated by the system at the time 
t can be defined as the fiux <I> of the Poynting vector S = (l/47r)E x B through 
a sphere S of radius R at the time tji, in the limit of large R. We have 



lim i?" 



df7„ S • n 



dflr, 

in 



[E X (n X E)] • n 



dQr. 

An 



(nxE)2, (8) 



where fin denotes the solid angle associated with the direction n, and integration 
is carried out over the total solid angle. Therefore the condition $ = is 
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equivalent to the requirement that the vector n x E vanishes for all directions 
n and times t: 



N 



= 



n X E = -n X > + v, 

^ (1 - n • \ 1 - n • Vi 



(9) 



The right-hand side of the above equation can be rewritten in a particularly 
compact form. In fact, using ^ we find that the partial derivative of U with 
respect to t at fixed n is 



dt 



„ l-v,(t,)-n 



We then have 



and 



9Zj 

dt 



1 - V,; • n 



(1 - n • V, 



1 - n • V,; 



(10) 



(11) 



(12) 



where all quantities with label i are evaluated at the time ti. Hence ^ can be 
rewritten as 



N 



nx 



= 0. 



Equation ([5]) implies that 



N 



N 



^ qi'Ziiti) = ^ qi{ti - t)n + C(n, t) , 



(13) 



(14) 



where C(n,t) is an arbitrary function such that C(n, t) • n = 0. Substituting 
the above expression into (fT5|) we then obtain 



C(n,t) = 0, 



so that 



C(n,i) =Co(n)+tCi(n), 



(15) 



(16) 



where Co(n) and Ci(n) are functions defined on the unit sphere |n| = 1, such 
that Co(n) • n ~ Ci(n) • n = 0. According to ([T]), we must have 



N 



N 

1=1 



(17) 



for all t. On the other hand, from ^1)1 and we obtain 



N 



> |Co(n)+iCi(n) 
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Therefore condition p?)) can be satisfied for i ^ oo only provided that Ci(n) — 
for all n. We then conclude that 

N 

"^q,[z,{U)~{U~t)n]=Co{n), (18) 

1=1 

with Co(n) • n = 0. We can express this result by saying that a bounded motion 
of a system of N charges does not radiate electromagnetic energy if and only if 
the quantity on the left-hand side of (fT5|) . where ti is determined by Eq. (O for 
alH = 1, . . . , iV, is independent of t for all unit vectors n. 

By differentiating with respect to t and using (fTU|) we obtain 

^ v.(t.)-[n.v(t)]n ^ (19) 

From the above formula it is easy to recover the well-known result that a single 
charged particle moving in a bounded region of space necessarily radiates. For 
= 1 in fact ([TO]) is equivalent to v(t) — [n • v(t)]n, which means that v(i) 
must be parallel to n. Since n can be varied independently of t, this condition 
implies v(t) = 0. Hence the particle must necessarily be static for all times in 
order to satisfy the condition of zero radiation. 



3 The system of two charges 

Let us consider the case N — 2. Condition of zero radiation can be written 
as 

V=(n-V)n, (20) 

where 

V = gi[l - n • V2(t2)]vi(ii) + (?2[1 - n • vi(ti)]v2(t2) . (21) 

In the above equation the times ti and ^2 are determined by Eq. ^ for i = 1,2. 
Hence we have 

t2-h = [z2{t2)-zi{ti)]-n, (22) 

which implies 

\t2 - h\ < \z2{t2) - Zi{h)\ . (23) 

In relativistic language, the above relation means that the two spacetime points 
(zi, ti) and (z2, t2), taken on the world-lines of particle 1 and 2 respectively, must 
have a space-like separation. Clearly, for any ti such that Z2(ti) 7^ zi(ti), there 
exists a finite interval of values of ^2, including the point t2 = ti, for which (j23p 
is satisfied. Note that, if zi{t) = Z2{t) for all t, then the two particles actually 
form a single compound particle with charge qi + q2, so that the situation is 
identical to the case of a single charge, which has already been considered at 
the end of the preceding section. Therefore, excluding this trivial case, in the 
following we shall always assume that zi{t) ^ Z2(t) for almost all t. 

Taking into account the arbitrariness of n and t, we see from (|22p that, if 
one takes any two times ti and ^2 satisfying ([25]) . then ([20]) must hold for all 
unit vectors n forming with Z2(i2) — zi(ii) the angle 

h - h 



|z2(i2)-Zi(ti)| 



(24) 
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The set of all such unit vectors forms a circle Cg of radius sin6' on the unit 
sphere. In order to satisfy Eq. 1^0]) . V has to be parallel to n for all n G Cg- 
However, for any 9 ^ 7r/2, the circle Cg is not contained in any plane containing 
the origin of the cartesian system. This implies, in particular, that almost all 
n € Cg do not lie in the plane containing vi(ii) and V2(i2)- Since V lies instead 
in this plane for any n G Cg, we see that Eq. ([20]) can be satisfied only if V = 0. 
From (PT|) it then follows that V2(i2) must be parallel to vi(ii). Substituting 
^2(^2) = Avi(ti) into the equation V = 0, and solving with respect to A, we 
obtain 

V2(t2) = (25) 

92 - [qi + q2)n ■ vi(ti) 

Let us initially suppose that qi + q2 — 0, which means that we are deal- 
ing with a neutral two-particle system (such as a hydrogen atom). Then 
becomes 

Vl(tl) =V2(t2) (26) 

for all t\ and t2 satisfying (|23|) . This means that, \it is such that Z2(t) ^ zi(i), 
then is satisfied for t\ — t and for all ti belonging to a finite interval 
containing t. This implies in particular that dv2(t2)/di2 — for = t. In 
the same way, by interchanging the role of particles 1 and 2, we also obtain 
that dvi(ti)/dti = for ti = t. We have thus proved that Vi(i) = V2(i) and 
Vi(t) = V2(t) = for all t such that Z2{t) 7^ zi(t). From this fact it easily 
follows that Vi{t) = V2{t) = constant for all t. Since the trajectories of the two 
particles were supposed to be bounded in space, we then conclude that 

vi(i) = V2(t) = foralH. (27) 

Therefore, for two charged particles such that qi + q2 = 0, the radiated power 
vanishes at all times only if the particles are static. 

Let us now consider the case qi + q2 7^ 0. If we take any two times ti and 
t2 satisfying ([23|) . we see from (f25|) and (|24p that n • vi(ti) must be a constant 
while n varies in the circle Cg. This implies that vi(ii) must be directed as 
22(^2) — zi(ti). By keeping t2 fixed and varying 0, one can actually prove that 
this fact is true for all ti belonging to a finite interval of time. A symmetrical 
result can also be proved for V2(i2)- We thus conclude that the whole motion of 
both particles must take place along a straight line. We have therefore reduced 
the problem to the study of a one-dimensional motion, and we will henceforth 
denote as 21 and Z2 the (scalar) coordinates of the two particles. We can rewrite 
(HI as 

[z2{t2)-zi{ti)]cose^t2~h, (28) 

and we have 

t2 — ti 

n • vi(ti) = vi{ti) cosO — vi{ti)- 



22(^2) - Zl{tl) ' 

with vi{ti) — dzi/dti, V2{t2) = dz2/dt2- Therefore from we obtain 
[Z2(t2) - ziihMqiViih) + (72^2(^2)] - (91 + 92)(i2 - h)viiti)v2it2) = 0. (29) 



This equation must hold for all ti and t2 satisfying Eq. ([23]) . By setting ti = 
t2 — t, where t is such that zi{t) ^ Z2{t), we obtain 

= qivi{t) + q2V2{t) = ^[qizi{t) + q2Z2{t)] , 
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which means that qizi{t) + q2Z2{t) is a constant. Since (?i + (72 7^ 0, by suitably 
choosing the origin of the z axis we can always set this constant to and obtain 

Z2(i) zi(t), V2{t) = -—vi{t) foralH. (30) 

92 92 

Hence, writing z and v in place of zi and vi respectively, Eq. (|29p becomes 

[912(^2) + q2z{ti)][v{t2) - v{h)] + (gi + q2){t2 - ti)w(ti)w(i2) = . (31) 
According to and ((501) . this equality must be true for all ti and t2 such that 
|92(i2 -ti)\< \qiz{t2) + 922(ii)| . (32) 
Interchanging ti and t2 in Eq- (PT|) . we also get 

[9l2(tl) + 92z(<2)]b(<2) - V{h)] + (91 + 92)(<2 - h)v{ti)v{t2) = (33) 

for 

192(^2 - tl)| < |9iz(ti) + q2z{t2)\ . (34) 

For all ti and t2 such that \t2 — ii| is sufficiently small, both conditions ((32l) 
and (jM)) are simultaneously satisfied. Therefore, subtracting ((55)) from (|5T|) we 
obtain 

(<72 - 9i)[^(t2) - 2(ti)]b(t2) - v{h)] = . 

If qi 7^ 92, the above equation implies that v(ti) = ^(^2) for all ti and t2 such 
that z(ii) 7^ z{t2)- But for such ii and t2 then (|3ip implies that v{ti) = v{t2) = 
0. On the other hand, for any regular function z(t), if z(ti) ^ z{t2) there must 
be a time i between ti and t2, such that z{i) ^ z{ti) and v{i) — z{t) 7^ 0. We 
see therefore that the hypothesis z{ti) ^ z{t2) leads to a contradiction. We 
must thus have z{ti) = z{t2) for all ti and t2, which means that the particles 
are motionless. 

Let us finally suppose that qi =92- Then (I5T|1 becomes 



[z{t2) + z{ti)][v{t2) - v{ti)] + 2{t2 - h)v{h)v{t2) - , (35) 

with z{t) = zi{t) = -Z2{t), v{t) = vi{t) = -V2{t). Let us divide dSS]) by - h 
and take the limit for t2 ^ ti = <. We obtain 

zit)z{t) + z^{t) = 0, 

or 

By integrating this equation we get 



z{t) ^y/a + bt, (36) 

where a and b are two integration constants. It is interesting to observe that 
this function z{t) is indeed a solution of Eq. ([55)1 for all a and b, so that it 
really describes a nonradiating motion of two equal charges. If ^ 0, by a 
suitable shift of the time axis it is always possible to set a = in ([36]) . Then 
the equations 

Zi{t) = Vbi, Z2{t) = -Vbi (37) 
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represent a nonradiating motion which is defined for all t such that bt > 0. 
If we also consider the conditions < 1, 1^2(^)1 < 1, which are imposed 

by special relativity, we must require that {l/2)y^\b/t\ < 1, or |t| > |6|/4. The 
motion is thus physically meaningful in the time interval 6/4 < t < +oo if 6 > 0, 
or— c»<t<6/4if6<0. However, in order for the motion to be bounded we 
must necessarily have 6 = in formula (j36p , so that we again conclude that the 
particles must be motionless. This result definitely excludes the possibility of a 
nonradiating bounded motion of any pair of arbitrary point charges. 

4 Discussion 

We have obtained a general condition for the absence of electromagnetic radia- 
tion from a system of moving point charges. This condition, which is expressed 
by formula p^ . seems at first sight very restrictive, since it must be satisfied 
for any arbitrary direction n. Using it we have deduced that two point particles 
of arbitrary electric charges cannot move for an infinitely long time within a 
finite region of space without radiating electromagnetic energy. However, an 
analogous result for more than two charged particles is at present not available. 
We have shown on the other hand that, if one only considers finite intervals 
of time (or, conversely, if one takes also unbounded trajectories into considera- 
tion), then nontrivial nonradiating motions of systems of point charges actually 
exist: Eq. ([57)1 provides an example of such a motion for qi = q2- 

In classical electrodynamics it is possible to remove in a relativistically co- 
variant way the divergences which are associated with the presence of point 
charges, and to obtain finite expressions for the energy and momentum of the 
complete system of particles and field [HI [17]. The conservation of these "renor- 
malized" quantities imposes on the particles the Lorentz-Dirac equation of mo- 
tion. However, the renormalized electromagnetic energy in the presence of point 
charges is no longer a positive definite functional of the field configuration [TB] . 
Therefore a spatially confined system of point charges might in principle keep 
radiating for an infinitely long time, while the electromagnetic energy contained 
in a finite volume including the particles diverges toward — cxi. For an isolated 
system, the electromagnetic field can be entirely expressed as a function of the 
dynamical variables of the particles at retarded times. Hence the divergence of 
the electromagnetic energy with increasing time must in any case be associated 
with an irreversible behavior of the system. It follows that a physically accept- 
able description of a stable system, such as an atom in its ground state, requires 
the existence of solutions which do not radiate, or which radiate at most a finite 
amount of energy during their whole history, starting from a given initial time. 

Let us consider a hypothetical confined solution of the Lorentz-Dirac equa- 
tion for two interacting particles with charges of equal modulus and opposite 
sign. Suppose also that this solution is such that the particles do not fall into 
each other either at finite or infinite times. Then the results of the preceding 
section suggest that, in order for the system to radiate at most a finite amount 
of energy, the accelerations of the two particles must tend asymptotically to 
zero. However such a motion is obviously not a solution of the Lorentz-Dirac 
equation, since the Coulomb attractive force does not asymptotically vanish. 
We conclude that the description of the hydrogen atom as an isolated system 
governed by the laws of classical electromagnetism is incapable of accounting 
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for the existence of bound noncoUapsing states. 

It is well known that, if radiation reaction is treated as a small perturbation 
of the mechanical trajectories for a charged particle in a Coulomb field, then a 
particle in a bound state should spiral toward the center of force and ultimately 
fall into it. At variance with the nonrelativistic case, the total energy radiated 
during such a process appears to be finite according to relativistic mechan- 
ics [IB]- The situation becomes however completely different if one treats the 
Coulomb problem by making use of the Lorentz-Dirac equation in an exact way. 
It has in fact been proved, either in the one-dimensional relativistic case [19] or 
in the three-dimensional nonrelativistic case [20], that there exists no solution 
of the Lorentz-Dirac equation for which the particle falls into the fixed center 
of force either at finite or infinite times. An analogous result, with not rigorous 
but quite convincing arguments, has also been obtained in the relativistic three- 
dimensional case [H]. It has also been shown in fTH] that, according to the 
relativistic Lorentz-Dirac equation, no collision can occur between two inter- 
acting particles of equal masses and opposite charges moving on a straight line. 
Let us now make the plausible hypothesis that these results can be extended to 
the case of two particles of different masses moving in three-dimensional space. 
In other words, let us suppose that for two particles there exist no collapsing 
solution at all. Since we have shown that nonradiating confined solutions of the 
Lorentz-Dirac equation do not exist, a confined solution should necessarily ra- 
diate an infinite amount of energy for infinite times. Hence the energy in a finite 
volume containing the system should diverge to — oo. Although we are unable 
at the moment to mathematically prove the impossibility of a noncoUapsing 
solution of this type, its existence would be quite surprising from a physical 
point of view. Taking all these facts into consideration, we are led to make 
the conjecture that, for the electromagnetic two-body problem with particles 
obeying to the Lorentz-Dirac equation, the only possible solutions are given by 
unbounded orbits. It is interesting in this respect to recall that, according to a 
recently obtained result [21j . only unbounded orbits can exist for a particle in 
a Coulomb field in the three-dimensional nonrelativistic case. 
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